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I. INTRODUCTION 



Random sequential adsorption (RSA) has attracted considerable interest due to its im- 
portance in many physical, chemical, and biological processes. In two dimensions (2D), RSA 
is a typical model for irreversible and sequential deposition of macromolecules at solid/liquid 
interfaces. Some examples of the wide range of applicability of this model include adhesion 
of colloidal particles, as well as adsorption of proteins to solid surfaces, with relaxation times 
much longer than the formation time of the deposit. A comprehensive survey on RSA and 
cooperative sequential adsorptions is given by Evans Recent surveys include Privman 
0,0], Senger et al. jj] and Talbot et al. j^. 

The simplest RSA model is defined by the following three rules: (i) objects are placed 
one after another in random position on the substrate; (ii) adsorbed objects do not overlap; 
and (iii) adsorbed objects are permanently fixed at their spatial positions (neither diffusion 
nor desorption from the surface are allowed). When the surface is saturated by adsorbed 
objects so that no further objects can be placed, the system reaches the jamming limit. The 
RSA models are broadly classified into continuum models and lattice models on the basis of 
the nature of the substrate. A quantity of central interest is the coverage 9{t), which is the 
fraction of the total substrate area occupied by the adsorbed objects at time t. Asymptotic 
approach of the coverage fraction 9{t) to its jamming limit 9 jam = 9{t —>■ oo) is known to be 
given by an algebraic time dependence for continuum systems 0, Q, [s], and by exponential 
time dependence for lattice models 0, [lo|, 11, 12|. 



An important issue in RSA is the influence of the shape of the depositing objects on 
the dynamics of irreversible deposition. RSA of many different geometric objects has been 
studied. For instance, Khandkar et al. [l3| have studied RSA of zero-area symmetric angled 



objects on a continuum substrate for the full range (0° — 180°) of values of the arm angle 
(f) and have observed that 9 jam — 9{t) ~ as expected. The value of the exponent a 
exhibits a crossover near = ° or 180 °, and is significantly lower in the case of the angled 
objects than in the case of needles. Wang and Pandey l4] have studied the kinetics and 
jamming coverage in RSA of self- avoiding walk chains on a square lattice. They reported 
that the growth of the coverage 9{t) to its jamming limit in the intermediate time regime 
can be described by a power law 9{t) ~ 9 jam — c/f . In contrast to theoretical predictions 
for the RSA of polydisperse object of regular shapes they find that effective exponent 7 
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depends on the chain length. They observed a crossover from a power-law variation of 
the coverage fraction 6{t) in the intermediate time regime to an exponential growth in the 



long time, especially for short chains. Budinski and Kozmidis [111 . Il2l | have carried out the 
extensive simulations of single-layer irreversible deposition using objects of different sizes 
and rotational symmetries on a square and triangular lattice. They fitted the asymptotic 
approach of the coverage fraction 6{t) to its jamming limit 6 jam by the exponential time 
dependence: 

eit)r^ejam-^0 exp(-t/a), (1.1) 

where and a are parameters that depend on the shape and orientational freedom of 
depositing objects [ll[ [3]. The shapes with the symmetry axis of a higher order have 
lower values of a, i.e., they approach their jamming limit more rapidly. The parameter A6 
decreases with the object size for the same type of shape. Later, by studying a reversible 



RSA of extended objects on a triangular lattice, Budinski et al. [15|] showed that the growth 
of the coverage 6(t) above the jamming limit to its steady-state value 6*00 is described by a 
pattern 6(t) = Boo — ^^-^/3[~(V''")^]) where Ep denotes the Mittag-LefHer function of order 
(3 G (0, 1). The parameter r is found to decay with the desorption probability P„ according 
to a power law r = A PZ^ ■ Exponent 7 is the same for all the shapes, but parameter A 
depends only on the order of symmetry axis of the shape. This confirms the crucial role of 
the geometrical character of the objects in deposition dynamics. 

In comparison to the irreversible deposition of pure depositing objects, very little at- 
tention has been given to the RSA of two or more species of different shap e and/or size 
although the latter problem is inherent in many experimental situations [igI . \\\ . The role 
of polydispersity has been studied numerically in a wide variety of conditions. Examples 
include, binary mixtures Il8^, \\^ , mixtures of particles obeying a uniform size distributions 

n n □ n 

|18l |. mixtures with Gaussian size distributions [18|, l20|, power-law size distributions |21| . 
etc. Theoretical works were restricted only to a binary mixtures of particles with very large 
size differences |22|. |23, 24 1. power-law size distributions {21! . 25|. or general continuous size 



distributions |26|]. The motivation of our present work comes from Barker and Grimson [27|, 
who investigated the adsorption of mixtures of lattice objects of different shapes but of the 
same size on a square lattice by means of computer simulation. 

In this paper we study the irreversible deposition of two-component mixtures of extended 
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objects on a 2D triangular lattice by Monte Carlo simulations. Simulations are performed 
for objects of various shapes. The depositing shapes are modeled by directed self-avoiding 
walks on 2D triangular lattice. A self-avoiding shape of length £ is a sequence of distinct 
vertices {ujq, . . . ,u!i) such that each vertex is a nearest neighbor of its predecessor, i.e., a 
walk of length i covers i+1 lattice sites. Examples of such walks for £ = 1, . . . , 6 are shown 
in Table [H On a triangular lattice objects with a symmetry axis of first, second, third, and 
sixth order can be formed. Rotational symmetry of order Ug, also called n-fold rotational 
symmetry, with respect to a particular axis perpendicular to the triangular lattice, means 
that rotation by an angle of 27r/ns does not change the object. Here we focus our interest 
on the influence of the order of symmetry axis of the shape on the kinetics of the deposition 
processes in two-component mixtures. 

The paper is organized as follows. Section [TTl describes the details of the simulations. 
We give the simulation results and discussions in Sec. IIIII Finally, Sec. IIVI contains some 
additional comments and final remarks. 



II. DEFINITION OF THE MODEL AND THE SIMULATION METHOD 



In order to make a systematic approach to this problem and to be able to compare the 
results for the mixtures to the results for the single-component systems, the simulations 
were also performed for pure shapes shown in Table I. The results are obtained for all the 
shapes that may have a different long-time behavior of 6{t) and different jamming coverages 
Ojam for the lengths of the walks i = 1, 2 and 3. The number of different shapes increases 



very fast with the length of the walk [28| and for the greater lengths we investigated a few 
characteristic objects for each length. We made at least two representative objects for each 
order of symmetry. 

At each Monte Carlo step a lattice site is selected at random. If the selected site is 
unoccupied, one of the six possible orientations is chosen at random and deposition of the 
object is tried in that direction. We fix the beginning of the walk that makes the shape at 
the selected site and search whether all successive i sites are unoccupied. If so, we occupy 
these i+1 sites and place the object. If the attempt fails, a new site is selected, and so on. 
After long enough time a jamming limit is reached when there is no more possibility for a 
deposition event. 
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It should be noted that the results in Table I differ from the estimates given in 12|]. In 
;he present work we deal with the above described conventional model of RSA, while in Ref. 

n 

12l | the end-on model of RSA has been used. In the end-on model the depositing object 
always checks all possible directions from the selected site. If the object cannot be placed 
in any of the six orientations, the site is denoted as inaccessible. The jamming limit for 
the end-on model is slightly larger than for the conventional model and the approach to the 
jamming limit is faster. 

The set of binary mixtures used in our simulations is shown in Table [TTl Each mixture 
(x) + (y) is composed of two lattice objects (A) - {U) from Table [H In the case of mixtures, 
at each Monte Carlo step a lattice site is selected at random, one of the objects that make 
the mixture is selected at random and deposition of the selected object is tried in one of 
the six orientations. If the attempt fails, a new site and a depositing object is selected at 
random. The jamming limit is reached when neither of the objects can be placed in any 
position on the lattice. 

The Monte Carlo simulations are performed on a triangular lattice of size L = 120. 
Periodic boundary conditions are used in all directions. The time is counted by the number 
of attempts to select a lattice site and scaled by the total number of lattice sites. The data 
are averaged over 500 independent runs for each shape and each mixture of shapes. 



III. RESULTS AND DISCUSSION 



For all the objects from Table I and for all the mixtures from Table II plots of 
\n{6jam — (^{t)) vs t are straight lines for the late stages of deposition. These results are 
in agreement with the exponential approach to the jamming limit of the form (1.1), with 
parameters a, A and 9jam that depend on the shape of the depositing object, i.e. on the 
combination of the objects making the mixture. 

Some of the plots for the pure objects are shown in Fig. 1. We can notice that there 
are four groups of lines with different slopes, corresponding to shapes with different order of 
symmetry. The values of the parameter a are determined from the slopes of the lines and 
they are given in Table I. According to a, all shapes can be divided into four groups: 

a) shapes with a symmetry axis of first order with a ~ 5.9; 
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b) shapes with a symmetry axis of second order with cr ~ 3.0; 

c) shapes with a symmetry axis of third order with a ^ 2.0; 

d) shapes with a symmetry axis of sixth order with a ^ 0.99. 

This means that the rapidity of the approach to the jamming hmit depends on the order of 
symmetry of the shape and the approach is slower for less symmetric shapes. The symmetry 
properties of the shapes have a crucial influence on the filling of small isolated targets on 
the lattice that are left for deposition in the late times of the process. Namely, a shape with 
symmetry axis of higher order has a greater number of possible orientations for deposition 
into an isolated location and an enhanced probability for the adsorption. 

Kinetics of the irreversible deposition of mixtures is illustrated in Fig. 2 where the plots 
of \n{9jam — 0{t)) vs t are given for some combinations of the shapes from Table I. These plots 
are straight lines for the late times of deposition, suggesting that in the case of mixtures the 
approach to the jamming limit is also exponential. In Fig. 2 lines with four different slopes 
can be noticed. The value of the parameter a is determined by the order of symmetry of 
the less symmetric object in the mixture and the values are: 

a) cr ^ 11.6 for the mixtures including an object with a symmetry axis of first order; 

b) cr ~ 5.65 for the mixtures that contain at least one object with symmetry axis of 
second order and no objects with symmetry axis of first order; 

c) cr ~ 3.97 for the mixtures that contain at least one object with symmetry axis of third 
order and no objects with symmetry axis of lower order; 

d) cr ~ 0.985 for the mixtures of the objects with symmetry axis of sixth order. 

In the late stages of deposition the less symmetric objects have to try all possible orientations, 
so they are responsible for the approach to the jamming limit. Deposition of mixtures is 
slower in comparison to the deposition of pure objects, except for the objects with symmetry 
axis of sixth order for which the parameter cr has the same values for the pure shapes and 
for the mixture. The retardation of the adsorption shows that each of the islands of the 
connected unoccupied sites is specific to a particular shape. If we imagine an island of 
unoccupied sites such that an adsorption of shape (x) is possible, prior to a succesful (x)- 
shape adsorption there will probably be some rejected attempts of (2/)-shape adsorption. 
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Therefore the overall process will proceed more slowly than each of the individual adsorption 
processes. 

Jamming configurations consist of clusters of objects of the same type. A few typical 
jamming configurations are shown in Figs. 3 a) - c). These clusters are more prominent 
for elongated and for symmetric objects. Namely, at very early times of the process the 
depositing objects do not "feel" the presence of the already deposited ones and are adsorbed 
randomly onto the surface. However, in the late stages of deposition the depositing objects 
must fit into small empty regions that favors the formation of clusters. If one examines the 
evolution of domains in the case of mixture (S) + (U) (Fig. 3 c)), it can be observed that 
growth of domains precipitated during early growth is substantially frustrated in the case 
of angled object {U). This is the cause of the severe limitation on the size of domains in the 
case of angled objects (U) as compared to those in the case of /c-mers (S). 

Value of Ojam for a mixture depends on the local geometry of the objects making the 
mixture. Values of the jamming coverage for different combinations of depositing objects 
are given in Table II. Qualitatively, we could say that it depends on the probability that the 
neighboring sites of an adsorbed object would be blocked by another adsorbing object from 
the mixture. For the same types of shapes the jamming coverage decreases when the size of 
the objects increases. 

For RSA of mixtures of fc-mers of various lengths ii and £2 on a square lattice it was 
found 29| that the jamming coverage for a mixture ofam^'^'' is always greater than either of 
the jamming coverages for the lengths ii and £2- Comparing the results from Table I and 
Table II we can see that for the mixtures of objects of various shapes this is not always the 
case. For a number of combinations of depositing objects given in Table II the jamming 
coverage for a mixture has greater values than the jamming coverages for the pure shapes 
making the mixture. However, there are also mixtures such as (E) + (F), (O) + (Q), (O) + 
(R), (P) + (Q), (S) + (T), (S) + (f/) that have a lower jamming coverage than one of 
the components. The jamming coverage of the mixtures is still greater than the jamming 
coverage of the other component. The mutual feature of these mixtures is that the jamming 
coverages of their components differ significantly. 

We have also performed extensive simulations in order to investigate the deposition pro- 
cesses for various fractional concentrations r^^^ and r^^^ of the shapes (x) and (y) in the 
reservoir, i.e. for various probabilities for choosing one of these shapes for a deposition at- 
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tempt. As an illustration, here we give the results for the jamming coverages and for the 
parameter a for a few combinations of the objects and for various fractional concentrations 
of depositing objects. In Figs. 4 a) - c) jamming coverages obtained for various composi- 
tions of mixtures (x) + [y) = (B) + (D), (S) + (T), (T) + {U) are shown vs the fractional 
concentration r^^^ (r^^^ = 1 — r^^^). We can see that 9 jam varies monotonously with the 
fractional concentration of one of the objects, growing with the concentration of the ob- 
ject with greater jamming coverage of the pure shapes. Depending on the combination of 
the objects and on their fractional concentrations it can be either greater than both jam- 
ming coverages of the pure shapes making the mixture or it can be lower than the higher 
single-component jamming coverage, but still higher than the other. When there is a large 
difference in the jamming coverages for the single-component depositions, such as in the 
case of the mixture (S) + (T) (Fig. 4 b)), the jamming coverage for the mixture can be in 
between these coverages, especially for the high fractional concentrations of the object with 
lower single-component 9 jam- 
Dependence of the parameter a on the fractional concentration r^^^ is shown in Fig. 5. 
The behavior of a differs from case to case. For example, for the mixture (B) + (D) (Fig. 
5 a)), a has a minimal value for fractional concentration r^^'^ ~ 0.6 and the process slows 
down when one of the fractional concentrations grows at the expense of the other. On the 
other hand, for the combination (T) + (f/), the value of a increases with the fractional 
concentration of the more symmetric object (T) (Fig. 5 c)). At first sight, this is somewhat 
surprising, since object (T) has a symmetry axis of sixth order and the lowest value of a 
in the case of a single-component deposition, and object (U) has a symmetry axis of first 
order and the highest value of a in the case of a single-component deposition. Indeed, the 
more symmetric objects reach their jamming coverage for a short time, but there are still left 
some small empty regions where the deposition of the less symmetric objects is possible. The 
lower is the fractional concentration of the latter objects, the longer is the time necessary 
to find all these empty regions left for the deposition. Despite these differences, the value of 
a has larger values for a mixture than for the pure shapes for all fractional concentrations 
of the components, so we can generally say that the deposition process is always slower for 
a mixture than for the pure shapes making the mixture. 
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IV. CONCLUDING REMARKS 



Kinetics of irreversible deposition of mixtures on a triangular lattice has been studied 
by Monte Carlo simulations. A systematic approach has been made by examining a wide 
variety of shapes and their combinations. The approach to the jamming limit was found to 
be exponential for all the shapes and all the mixtures. The rapidity of the approach depends 
only on the symmetry properties of depositing objects. For two-component mixtures the 
value of the parameter a is determined by the order of symmetry of the less symmetric 
object in the mixture. 

It was shown that the kinetics of irreversible deposition fastens with the increase of the 
order of symmetry of the shape. On the contrary, the adsorption- desorption pro cess of 
asymmetric shapes is faster than the same process of more symmetric shapes [15|. The 
reason for both effects lies in the filling of small empty regions in the late stages of the 
processes. In the case of irreversible deposition, objects with symmetry axis of higher order 
have greater number of possibilities for deposition into fixed small empty regions in the 
late stages of deposition, while in the reversible case the enhanced rate of single particle 
adsorption prolongates the approach to the steady state coverage. Namely, in the reversible 
case, when the value 9 jam is reached, the rare desorption events are generally followed by 
immediate readsorption. These single-particle events do not change the total number of 
particles. However, when one badly placed particle desorbs and two particles adsorb in the 
opened location, the number of particles is increased by one. These collective events are 
responsible for the density growth above 9 jam- The increase of the order of symmetry of the 
depositing object enhances the rate of single-particle readsorption. This extends the mean 
waiting time between consecutive two-particle events and the approach to the steady state 
coverage is slower. 

Jamming configurations for RSA of mixtures consist of clusters of blocked sites and 
of clusters of objects of the same type. The dimensions of these clusters are greater for 
elongated and for symmetric objects. Jamming coverage for a mixture is always greater 
than the jamming coverage of the component with lower 9 jam and it is often greater than 
either of the jamming coverages of the components making the mixture. 

When the fractional concentration of one component grows at the expense of the other 
component in the mixture, the jamming coverage varies monotonously. It increases when 
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the concentration of the object with greater single-component jamming coverage increases. 
On the other hand, the rapidity of the approach to the jamming coverage doesn't always 
show a monotonic behavior. For many combinations there is a minimal value of a, i.e. the 
jamming limit is reached in the shortest time, for a certain value of fractional concentration 
of one of the components making mixture. 
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Figure 1 

(Color online) Plots of ln{6jam — ^(^)) vs t for various objects from Table I. Different slopes 
correspond to different symmetry orders. 
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Figure 2 



(Color online) Plots of \n{6jam — vs t for various mixtures from Table II. The slopes 
are determined by the order of symmetry of the less symmetric object in the mixture. 
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Figure 3 a) 



(Color online) Snapshots of patterns formed during the RSA of mixture {B) + {D) correspond 
amming state; (5)-red, (D)-green. 
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Figure 3 b) 



(Color online) Snapshots of patterns formed during the RSA of mixture {S) + {T) correspond 
amming state; (S')-red, (T)-green. 
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Figure 3 c) 



(Color online) Snapshots of patterns formed during the RSA of mixture (S) + (U) correspond 
amming state; (S')-red, (?7)-green. 
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Figure 4 a) 

Jamming coverage for the mixture {B) + (D) vs the fractional concentration r^^\ The dotted 
and the dashed hues represent the jamming coverages for the single-component deposition 
of shapes {B) and {D), respectively. 
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Figure 4 b) 

Jamming coverage for the mixture {S) + (T) vs the fractional concentration r^^\ The dotted 
and the dashed hues represent the jamming coverages for the single-component deposition 
of shapes (5) and (T), respectively. 
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Figure 4 c) 

Jamming coverage for the mixture (T) + (U) vs the fractional concentration r^'^\ The dotted 
and the dashed hues represent the jamming coverages for the single-component deposition 
of shapes (T) and (U), respectively. 
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Figure 5 a) 

Values of the parameter a for the mixture {B) + {D) vs the fractional concentration r'^^\ The 
dotted and the dashed lines represent the values of a for the single-component deposition of 
shapes {B) and (D), respectively. 
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Figure 5 b) 



Values of the parameter a for the mixture {S) + (T) vs the fractional concentration r^^\ The 
dotted and the dashed lines represent the values of a for the single-component deposition of 
shapes {S) and (T), respectively. 



22 







- 


■ mixture (T)+(U) 


■ 




pure (T) 






- - ■ pure (U) 




- 
- 


■ 






■ 




■ 


■ 




1 1 1 1 1 1 1 1 1 1 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 5 c) 



Values of the parameter a for the mixture (T) + {U) vs the fractional concentration r^^\ The 
dotted and the dashed lines represent the values of a for the single-component deposition of 
shapes (T) and (?7), respectively. 
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(Q) 

(R) 



(S) 
(T) 
(U) 



A 



y 



A 



A7 
\ 



6 



jam 

0.9139 



0.8362 
0.8345 
0.7970 



0.7886 
0.7653 
0.7739 
0.7404 

0.7651 

0.7226 
0.7593 

0.7744 

0.7742 



0.7587 



0.7372 
0.7210 
0.5740 
0.6758 



0.7212 
0.6695 
0.6443 



A9 



a 



0.105 3.12 



0.134 2.94 
0.0813 5.78 
0.115 1.96 



0.0944 3.12 

0.0813 5.74 

0.0602 5.98 

0.0973 3.03 

0.0805 5.75 

0.0573 5.84 

0.0590 2.93 

0.0578 6.04 

0.0608 5.96 



0.115 2.87 



0.0837 3.02 
0.0944 1.99 
0.0707 0.986 
0.0829 5.81 



0.0954 2.84 
0.0773 0.994 
0.0721 6.07 



TABLE I: Parameters A6 and a determined using Eq. (jl.ip for various shapes (x) of length ^ on a 
triangular lattice. The colors (online only) are associated with the different order ni'^^ of symmetry 
axis. The typical statistical errors are estimated to the last given digits. 
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shapes 



(^•) , {y) 



jam 



{A 
{A 
{A 
{A 
{A 



+ {B) 
+ {E) 

+ {0) 
+ {S) 



+ 



+ . 



+ . 



2 + 2 

2 + 2 

2 + 2 

2 + 2 

2 + 2 



1 + 2 

1 + 3 

1 + 4 

1 + 5 

1 + 6 



0.9202 
0.9191 
0.9196 
0.9195 
0.9198 



0.5401 0.3801 
0.5833 0.3358 
0.6144 0.3052 
0.6143 0.3052 
0.6553 0.2645 



0.125 
0.129 
0.131 
0.135 
0.157 



5.52 
5.66 
5.77 
5.73 
5.76 



{C 



+ {C) 
+ {D) 
+ {D) 

+ {P) 



+ 

.+ A 



^+ A 
A + 



2 + 1 

2 + 3 
1 + 3 

3 + 3 



2 + 2 

2 + 2 

2 + 2 

2 + 5 



0.8526 
0.8591 
0.8624 

0.8211 



0.4191 0.4335 
0.4330 0.4261 
0.4406 0.4218 

0.4946 0.3265 



0.0365 
0.0781 
0.0493 

0.0983 



11.42 

5.51 

11.62 

3.95 



+ {F) 
+ {K) 
+ {K) 



+ 



V 



_+ A7 

y+ A7 



2 + 1 
2 + 2 
1 + 2 



3 + 3 
3 + 3 
3 + 3 



0.7876 
0.8109 
0.8140 



0.3903 0.3973 
0.3917 0.4192 
0.3960 0.4180 



0.0354 
0.0620 
0.0523 



11.84 
5.87 
11.55 



(O 

(O 
{P 
{P 
{Q 

(Q 



+ ip) 

+ {Q) 
+ {R) 
+ {Q) 
+ {R) 
+ {R) 

+ {T) 



+ 



+ 



2 + 3 



2 + 6 



2 + 1 



3 + 6 



3 + 1 



6 + 1 



6 + 6 



5 + 5 



5 + 5 



5 + 5 



5 + 5 



5 + 5 



5 + 5 



5 + 6 



0.7647 



0.6791 



0.7086 



0.6833 



0.7485 



0.6822 



0.6222 



0.3415 0.4232 



0.3608 0.3183 



0.3612 0.3474 



0.3998 0.2835 



0.4156 0.3329 



0.3227 0.3595 



0.2873 0.3349 



0.0639 



0.0534 



0.0460 



0.0686 



0.0420 



0.0429 



0.0758 



5.41 



5.54 



11.21 



3.98 



11.74 



12.07 



0.985 



(S 
(S 
(T 



+ iT) 
+ {U) 
+ {U) 



,+t:7 



+ 



+ 



2 + 6 
2 + 1 
6 + 1 



6 + 6 



6 + 6 



6 + 6 



0.7125 



0.6833 



0.7087 



0.3218 0.3907 



0.3533 0.3300 



0.3955 0.3132 



0.0344 



0.0437 



0.0450 



5.76 



11.22 



11.57 



TABLE II: Coverage fraction ofj^^^^ for various binary mixtures (x) + (y) of shapes A - V from 
Table I. The colors (online only) are associated with the different order ni^^ of symmetry axis. The 
typical statistical errors are estimated to the last given digits. 
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